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Abstract. 

We investigate aspects of axion as a coherently oscillating massive classical scalar 
field by analyzing third order perturbations in Einstein’s gravity in the axion-comoving 
gauge. The axion fluid has its characteristic pressure term leading to an axion Jeans 
scale which is cosmologically negligible for a canonical axion mass. Our classically 
derived axion pressure term in Einstein’s gravity is identical to the one derived in the 
non-relativistic quantum mechanical context in the literature. We show that except for 
the axion pressure term, the axion fluid equations are exactly the same as the general 
relativistic continuity and Euler equations of a zero-pressure fluid up to third order 
perturbation. The general relativistic density and velocity perturbations of the CDM 
in the CDM-comoving gauge are exactly the same as the Newtonian perturbations to 
the second order (in all scales), and the pure general relativistic corrections appearing 
from the third order are numerically negligible (in all scales as well) in the current 
paradigm of concordance cosmology. Therefore, here we prove that, in the super-Jeans 
scale, the classical axion can be handled as the Newtonian CDM fluid up to third order 
perturbation. We also show that the axion fluid supports the vector-type (rotational) 
perturbation from the third order. Our analysis includes the cosmological constant. 
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The currently popular concordance cosmology based on perturbed Friedmann world 
model demands the presence of substantial amount of non-baryonic dark matter. Despite 
some problems concerning small-scale (less than some kpc) clustering properties, the cold 
dark matter (CDM) still plays a popular role as the dark matter. Axion, notwithstanding 
its still hypothetical nature, is known to be one of the prime candidates for the CDM 

[I]- 

Axion as a coherently oscillating scalar field is known to behave as a pressureless 
fluid in the background world model [1], and to the linear order perturbation in the 
Newtonian context [2]. The CDM nature of linear order perturbation of the axion in 
Einstein’s gravity was studied in [3l HJ [5l |6l [7] . The CDM nature of axion to the linear 
order was shown in various gauge conditions: these are the zero-shear gauge mn, the 
synchronous gauge [1], the uniform-curvature gauge [5], and the axion-comoving gauge 
[7]. Only in the axion-comoving gauge the axion is shown to behave as the Newtonian 
CDM in all cosmological scales [7]. The axion, however, has its characteristic pressure 
term with Solar-System size Jeans scale for the canonical axion mass [21 El Eld, see Eq. 
(jHTll . We may increase the Jeans scale to cosmologically relevant one by reducing the 
axion mass; roles of this extremely low mass axion as a potential variant of warm dark 
matter with small scale cut off in the density power spectrum deserve further study [8]. 

Previously we have shown that the axion as a classical massive scalar held behaves 
as the CDM to the second-order perturbations in all cosmological scales [9]. Except 
for the axion pressure term relevant in the small-scale limit, the relativistic continuity 
and Euler equations of the axion huid are shown to be exactly the same as Newtonian 
equations to the second order in all cosmological scales including the super-horizon scale; 
in the case of a zero-pressure and irrotational huid, see [TO] . 

Here, we extend our study to the third-order perturbation in the axion-comoving 
gauge. The leading nonlinear contribution to the density and velocity power spectra 
needs the third order perturbation. We will consider the axion as a classical massive 
scalar held in Einstein’s gravity. We will use the fully nonlinear perturbation formulation 
in mm We will show that the general relativistic CDM correspondence continues 
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to the third order with due presence of the axion pressure term in the small-scale: for 
the relativistic continuity and Euler equations of axion, see Eqs. fH6l) and Wi\ . We 
show that the axion pressure term classically derived in Einstein’s gravity (to third 
order) coincides exactly with the one derived from non-relativistic quantum mechanics 
in the literature. We show that from the third order the axion generates the rotational 
perturbation: see Eq. fl57|l . For extremely low mass axion where the axion pressure term 
becomes cosmologically important, we note that the third order pure general relativistic 
correction terms are also affected by the axion pressure term: see Eqs. (155)1 and (|56il . 

Readers who want to skip the technical details of deriving the general relativistic 
continuity and Euler equation for axion fluid based on relativistic nonlinear perturbation 
analysis can go directly to Secs. |5] and [6] after reading Erst two paragraphs of Sec. 01 
For the beneht of such readers, below we briefly explain the basic method of analysis in 
the case of axion where we need fluid-held mixed formulation of nonlinear cosmological 
perturbation with proper time averaging. 

As the axion we consider a classical minimally coupled massive scalar held and take 
temporal averaging over the coherent oscillation. The Einstein’s equation becomes 

Ga, = 87rG{f,b), ( 1 ) 

where we take the temporal averaging for the scalar held to get the averaged energy- 
momentum tensor, and the averaged huid (density, pressure and velocity) quantities; 
although the average should be applied to the held quantities, here only to show the 
method we impose it to the energy-momentum tensor. For the energy conservation 
equation we can use 

(f„y = 0, or = (2) 

where the latter expression is the ADM energy-conservation with ha the normal-frame 
four-vector. We note that in the case of axion it is troublesome to use the covariant 
energy conservation = 0 as the huid four-vector also involves the scalar held 

[HIE]. For the massive scalar held Tab involves the scalar held in quadratic combination 
only, see Eq. ([8]). For the equation of motion we do not take average, see Eq. flTSl) . As 
the axion-comoving gauge we impose the scalar-type (longitudinal) perturbation part 
of (T°) equals to zero to all orders in perturbation, see Eq. (El]). We will use the fully 
nonlinear equations expressed using the huid quantities mm; one important point in 
the analysis of axion is to hrst express the huid quantities in terms of the scalar held, 
and then take average over the quadratic combination of scalar held. 

We set c = 1 = h except for Sec. [5l 

2. Metric and fluid quantities 

We consider scalar- and vector-type perturbations in the Friedmann world model with 
the general background curvature, K. Our metric convention valid to the nonlinear 
order is [13] 

{1 + 2a) drf — 2a^{l3^i + B^^^)dr]dx^ 
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(1 + 2^p) 'jij + 2'y^i\j + + C^^|/ + 2CO dx'‘dx\ 


-<0) 


N) 

■'ij 


(3) 


where a is the cosmic scale factor; a, /9, 7 , (^, Bf\ and C\f are arbitrary functions 
of space and time with = 0 = for the vector-type perturbations and 


= 0 = - for the tensor-type perturbation; the spatial indices of perturbation 

variables are raised and lowered by 7 ^^ as the metric, and the vertical bar indicates 
the covariant derivative based on 7 jj as the metric; for concrete forms of 7 *^, see [ 12 ]. 
The scalar-, vector- and tensor-types decomposition is generally valid considering fully 
nonlinear perturbations [T3|. To the nonlinear order the three-types of perturbations 
couple with each other in the equation level. In this work we ignore the tensor-type 
perturbation . 

The energy-momentum tensor in the energy-frame = 0) is given as 

Tab S'UaUb T P (^Qab T UqU})^ S T^abi (4) 


where Jl, p, Ua and itab are the energy density, pressure, fluid four-vector, and the 
anisotropic stress, respectively [I5]; we have u'^Uc = — 1 , ^ah = ^ba and = 0 = 

We introduce the perturbed order fluid three-velocity Vi as m 

Ui = a^Vi, 7 = ^ . (5) 

I _ jihk. 

V 1 + 2 V’ 

The index of Vi is raised and lowered by as the metric; we ignore the anisotropic 
stress TTab as the minimally coupled scalar held does not support the anisotropic stress. 
As the spatial gauge condition we take 

7 = 0 = Cf>, (6) 

on the metric tensor. This is the only spatial gauge condition which (together with 
our temporal gauge condition) allows the remaining variables to be gauge invariant 
to the fully nonlinear order [131 IIHl HI]- With this spatial gauge condition we set 
Xi = a{/3^i -|- = x,i + Xi^\ and the metric becomes 

(1 -I- 2a) dp^ — 2aXidpdx^ -|- (1 -|- 2(^) Xijdx^dxT (7) 

This allows us to expand the cosmological perturbations to fully nonlinear order in an 
exact manner mm- 


3. Minimally Coupled Scalar Field 


A minimally coupled scalar held is given as 


Tab — 


Pa(P,b 




9ab- 


In the energy-frame we can show TTab = 0 mm- 
The huid quantities follow from Eqs. (Hj) and 


p = ^(j) +V, p = ^0 - F, 


Vi = 


( 8 ) 


[see Eq. (5.6) of 


070 


(9) 
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where 


— 0.rX^ — 


= — — + 


X 


1 D(j) Dcf) 

^ Dt’’ Dt J\f \dt ' a^(l + 2lp) 

1 


V,; 


7 = 


1 - 




Af = 


\ 


1 + 2a + 


x’^Xk 


0^(1 + 2(p) 


a^(l+2<p){D(t>/Dt)^ 

The equation of motion, 0’^^^ = dV/{d^), gives [see Eq. (5.10) of 


+ (3HAf - AAk - ^ - 






Af a^Afi^l + 2(p) _ 




x'x^ 


(j) + Af‘ 


4’,i\j + 


dV 


+ 


2 x* 


d(j) a^(l + 2lp) 

Af^ (Af\^ 


a^{l + 2(f) y a?{l + 2(f) 

(‘^TjKr-Kr ^ ^ 

-(^3 AA Afn ^ aW(l + 2(^)j a2(l + 2v?) 


+ 


ipy 


a?{l + 2^)\Af 1 + 2^ 


X 


X 


X 


a2(l + 2ip)) 0^(1 + 2if) V1 + 2ip^ 


|fcJ 


= 0 . 


( 10 ) 


( 11 ) 


4. Axion 


We consider the axion as a massive scalar held with a potential V{(f)) = ^mScfA. The 
pseudo nature of axion is not relevant in our cosmological consideration. We consider 
only the classical nature of axion; we are not capable of participating in the ongoing 
controversies about the potential roles of quantum nature of the axion as a Bose-Einstein 
condensation 0[16]; see below Eq. (15^ though. We have 


H 

m 


2.133 X 10"2®h 


/ lO-^eV \ 
\ m j 



( 12 ) 


where Hq = 100hkmsec“^Mpc“^ is the present Hubble parameter with H = a/a. As 
the axion coherently oscillates we strictly ignore H/m higher order terms. 

We decompose the held and huid quantities to the background and perturbed parts 


as 


0 = 0 + 50, /i =/i + 5/i, p = p + 5p. (13) 

We have a background solution mu 

(j){t) = [0+0 sin {mt) + 0_o cos {mt)] , (14) 

where 0+o and 0_o are the constant coefficients. We take average over time scale of 
order m~^ for all huid quantities associated with the axion as a oscillating scalar held. 
|1] . We have [H |5] 

+ = (0^0 + <P-o) ) P = 0- (15) 

Thus, for the background the axion evolves exactly the same as a pressureless ideal 
huid pp. This conclusion is valid in the presence of both the spatial curvature and the 
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cosmological constant A in the backgronnd, and in the presence of other flnids. For the 
nonlinear pertnrbation, here we consider a single presence of the axion field; considering 
the presence of other flnids and fields is trivial and those are not affected by the special 
natnre of axion flnid. 


4-1. Fluid quantities and equation of motion 

Inclnding pertnrbations we expand mn 

f = (j) + 6(j) = + <50+) sin (mt) + + (50_) cos (mt), (16) 

to fnlly nonlinear order, where <50+ and (50_ are arbitrary fnnctions of space and time. 
In the analysis we ignore H/m higher order terms with 50+ ~ H6(j)±. 

For the flnid qnantities, from Eq. ([9]), to the fnlly nonlinear order, we can show 




P 




(0 ) + 2 


A 


0^(1 + 2ip) 


(00,i) + 




a\l + 2^y 




-± -m^( 0 ^) 

2a^{l + 2^y^ 2 ' 


= m“- I —^ ± 1 
4 VAA2 


,-3 


+m 


X 


+ 0—o) A 2<^ (0+o<50+ + 0—o<50—) + <50+ + 50?_ 

(0+o50-,i - 0_o50+,i) + 50+50-,* - 50_50+, 


2 aW 2 (l + 2y L 

(50+,*50+,j + 50_,*50_,j) 


A*A-^ 


4aW2(l + 2ipy 

~ 4+(l + 2y) {‘’S'S't+.i + , (17) 

where the npper and lower signs correspond to the J1 and p parts, respectively. In the 
temporal comoving gange (hypersnrface) condition we will take, the n* part is hxed by 
the gange condition, and later we will identify other variable as the velocity pertnrbation, 
see Eqs. (|3^ and (03]). In the axion case where we need proper averaging, the parameter 
Vi becomes ambignons, see paragraph containing Eq. (ED. 

The eqnation of motion in Eq. flTTll . to the fnlly nonlinear order, gives 

(TV — 1) =F3Ffm0:po + (1 +-A/") 0+0 

=F2m50+ =F 3ifmA/'50+ + 50+ + 3HJ\f6f± + 50+ 


Af K T —p T 




-3/2 


2 x* 


0^(1 + 2(p) 

Af^ 


Af a‘^Af{l + 2(p) _ 


(t"i50zp,* + 50+,*) 


=F"i 0 +o 
1 , 


if 0+0 =F m50+ + 50d 


0^(1 + 2(p) 






0^(1 + 2ip) 


^ 0 ±,i|j 


+ 




+ 


C' 


'{l + 2p) {AT l + 2ip 


























CONTENTS 


7 




X 


[^HU Uk ^ a?U{l + 2^)) a?{l + 2^) 


X 


+ 


X 


X 


^ 0 ±,i 


X 


sin {mt) 
cos {mt) 


= 0 , 


( 18 ) 


^( 3^(1 + 2(f)j a^{\ + 2(p) \ 1 + 2(pj 1 ^ 

where the upper and lower signs correspond to the sin and cos parts, respectively; here 
we have not ignored i7/m-order terms yet. 

4 . 2 . Axion-comoving gauge 
We can decompose 

Vi = -v^i + vf\ = 0. (19) 

In the fluid case the comoving gauge takes 

7 = 0 , ( 20 ) 

as the temporal hypersurface (slicing) condition. For a scalar held the comoving gauge 
is the same as 5(f) = 0, the uniform-field gauge to fully nonlinear order mm- 

In the axion case, however, as we have to take proper average over the oscillating 
held, the concept of Vi and its decomposition become ambiguous. Thus, for the temporal 
gauge condition we go back to the original energy-momentum tensor and take the scalar 
part of T? equal to zero as the comoving gauge. Thus, in the axion-comoving gauge we 
can impose 

(?“)'• = 0 , ( 21 ) 

to all perturbation orders. Together with the spatial gauge condition in Eq. (|6]) the 
axion-comoving slicing condition completely Axes the gauge mode. Each perturbation 
variable in these conditions has a unique gauge-invariant counterpart to all orders in 
perturbation [T3 l ITOt [TT|. 

Using the inverse metric presented in [12], from Eq. 


^^ ^ ^ we can show 

1 

{(f>+o5(f>-,i - (f>-o5(f>+,i) + m {6(f>+6(f>-^i - 6(f>-6(f>+^i) 

— - {5(j)+^i5(j)+j + 5(j)-^i5(j)-j) 


1 


2aAf^ 


( 22 ) 


0^(1 2 ( p ) 

As we mentioned below Eq. ([T|) we can omit average notation for Tab- Up to this point 
all equations are valid to fully nonlinear order. 

Now, we consider perturbations to the third order. From the gauge condition in 
Eq. we can show 

6(p+ 6(p- a _3f2^lo + <P-o ^-lv 7 i fxj. xj. A 


'+ 0«'-0 


(23) 
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to the third order, Thus 

=-^ \6ct>+,6ct>+,iX^ - V.A-V^' (<50+,,50+,fcx' 

m (pj^Q 

which is pure third-order vector-type perturbation. In the fluid formulation Eqs. 
and dS]) give 


(24) 


T)' = (fi + p) 

By introducing Uj in this fluid relation and using Eq. flT^ . we have 

u = 0, 


(26) 




_ _ rj-^y _ 




+0 


5ct>+,5ct>+,jX^ - V.A-V^' . (26) 


Thus, axion has the rotational perturbation appearing from the third order; from Eq. 
dS]) we have Ui = avf\ 

To the third order in the axion comoving gauge Eq. (Ell can be arranged as 


T 

P 


= T- 


l±^f^ 

2Af^ 


,3/2 


1 + 2 


0+0 + 0-0 
3 


(0+0(50+ + 0-0(50-) + -7^(50: 


000 


These are the energy density and pressure supported by the axion fluid to the third- 
order perturbation. Besides the Einstein’s equations to be presented below, we need an 
expression of SJ\f in terms of the axion field. This relation follows from the equation of 
motion, and will be presented in Eq. fl35l) . 


4-3. Einstein equations 


Equations valid to fully nonlinear perturbation orders are presented in [HI |T2] . Here 
we will again borrow some equations from these works. 

As explained in the paragraph containing Eq. dlD we can use the equations expressed 
in terms of fluid quantities but have to change fluid quantities to the scalar held before 
we take average over the oscillating held. After averaging, we will strictly ignore H/m- 
order terms, and will keep only first order in A/{m?a?). It is convenient to have 

K = S H6, Xi = X,i = (28) 

to the linear order. 

From the covariant energy conservation and the trace of ADM (Arnowitt-Deser- 
Misner) propagation equations in Eqs. (3.8) and (3.4) of [12], respectively, to the third 
order, we have 


/i + 


0^(1 + 2(p) 


Sp,kX^ + A/" (p + p) {3H - k) = 0, 


(29) 


-3 


AA' ' + + 


47tG 


A' 

(h + 3p) - - 


1 . „(i AA/" 

+ "77^ + 2 — K H-- , ——-—r 

AA a a^M{l + 2(p) 
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1 + 2 (^ 


+ T^,Kl 


(30) 


where from Eq. (3.10) in 


K K- = - , 

' * aW2(l + 2(y9)2 1 2 


{Xi\j + Xj\i) -\x"\iX^y 


^XV'^ {Xi\j + Xj\) - ^X*|*xVjj 
to the third order. A perturbation variable k is dehned as 


Kl = -3- + K = = -hh 


;c) 


(31) 


(32) 


where Kij is the extrinsic curvature, and 0*-”^ is the expansion scalar of the normal- 
frame four-vector; later we will identify k as the divergence of velocity perturbation, see 
Eq. fj43l) . We also need the ADM energy constraint and ADM momentum constraint 
equations to the second order. From Eqs. (3.2) and (3.3) in [12], respectively, we have 

= -A-kGqS - -k + ^ (AipNp + - 3(/9 (1 - 2p) ^ 


a 






Aa‘^ 


\x"^^ (Mb' + Xi|i) - 

^ (Am + - \x'"\ki 


2 

Q^,i + 


aW(l + 2(p) [2 




Ik 


+ V (x'V,i + xt^‘ - -^six v.k 


(33) 


(34) 


valid to the second order. 

4 . 4 . Results 

Now we assume a flat background curvature, K = 0. From the equation of motion in 
Eq. (ITSll we can show 


hAA = (1 - <h + 


A 


where $ = 50+/(a“^/^0+o), and we set Af = l + 6Af. This relation follows from dividing 
Eq. ffTSj) by 0=^0 and summing the two relations. 

From Eq. 0271) we have 

3 = 24+ 4=-(1 + 4) ^4-^4%. 

— = - (1 + 4) 


* “ + liv) (■’" + S'’”) + idw 


1 - I S-'6., 


(36) 


Thus, we have 

6Af= ^ 


Am? 


f 1 c- 

3 ,n\ A , 1 

/ 1 A 

A . 1 A ,0 

1 - -5 + 

-(5M -5-- 

1 

1 

—5^ + — -5^ 

V 2 

^ a? 4 

\ 2 J 

a? 8 


1 AV 1 + ^ 
2 m 2 a 2 x/r^T’ 
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5p 

fi 


Am? 


_ ( 1 + _ I 52 

2 8 


1 

^5 + T 






8 a2 


2a? 




,(37) 


where the compact expression in dM is motivated by the known expression in the 
literatnre on non-relativistic Bose-Einstein condensation in Minkowski space-time ra. 
see below Eq. (|52]). The ADM momentnm conservation in Eq. (3.7) of [12] gives 


p,i + [t + p)^ = -T 


M ^m‘^0? 




(38) 


and this is consistent with Eqs. fl3^ and 

Now, Eqs. (Enj), flSU]) and flMD . respectively, give 


K 


— (1 - 2(p) 5 iX* - 5 k = 0, 


(39) 


k + 2Hn — A7iGfi6 -I- ^ (1 — 2lp) k^jX* — 


(1 - 4(p) 


iXij + Xj,i) - o i^xY 


+ 


X"T’^X,ij 


;X"tAx 


A 


5M = 


A A^/m 

2m?a‘^ VTTT^ 


valid to the third order, and 


/ A \ O Ll 

+ ~x). + 7 — 


3 A (t,) _ ■ 


A a? 


2 

(2(pAx - T,jX^) . + - {T,ijX^ + xAt) 


(40) 


(41) 


valid to the second order. We note that compared with the CDM case in [18] the only 
correction from the axion natnre appears in the right-hand-side of Eq. (HOjl . 

From Eq. (1TT|) we have xY'^ = 0 to the linear order. To the second order Eq. (ITT|) 
can be decomposed as 


A 1 
«+ ^ 




2ipAx - (p,iX' + V* {(p,ijX^ + xAt) 
T,ijX^ + xAt - ViA" {pdkX^ + Xj^t) 


^ (v) _ 

9 Xi 9 

A 1 1 

~Xi = + 7 ^i- 


a 


a 


a 


a 


(42) 


5. Relativistic continuity and Euler equations for Axion 


Now we have complete relations for axion flnid valid to the third order pertnrbation in 
Einstein’s gravity. We will arrange the eqnations so that we can readily compare with 
the Newtonian flnid eqnations. 

We identify 6 as the density contrast 6g/g, and k a pertnrbed part of trace of 
extrinsic curvature in Eq. fl52]) as the divergence of velocity perturbation [THl US] 


K =-V ■ U =- U. 


(43) 


a a 

Here we assume the flat background curvature but include the cosmological constant. 
From Eq. fl42p . to the second order, we have 


-X = u + A-'(VX + Y), 
a 


(44) 
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X = 2ipV ■ u — u ■ V(p + -A ■ (u ■ VV(^ + uA(^), 

Y = 2 u ■ VV 93 + uA^? — VA~^V ■ (u ■ VV</9 + uA(y9) , (45) 

with V ■ Y = 0; X* cind (f are metric variables in Eq. ([7]), interpreted as the perturbed 
shear (of the normal four-vector flow) and perturbed spatial curvature, respectively 

III1II2]. 

Equations fl3^ and fl40|) give the general relativistic version of the continuity 
(mass-conservation) and Euler (momentum-conservation) equations valid to third order 
perturbations in the axion case. These are 


5 + -V ■ u + -V ■ (5u) = - (V5) • \2Lpvi - A’^ (VA + Y) 
a a a ^ 


-V ■ (u -I- -u) -I- 4 :7iGg6 H —rV ■ (u ■ Vu) — 

a \ n / 


h^A A^/TTS 
2mSa‘^ + 5 


= - -(p\i ■ V (V ■ u) -h 4V • 




U ■ Vu — gUV ■ u 


(46) 


-.YV ■ u + u ■ (VV + Y) - A 

3 


u-A-^(VX + Y) 


In order to close the equations we need (p valid to the linear order. From Eq. 
have 


.A 


d 1 , 


(47) 


we 


(48) 


-ip = —AnCgd -\ -V ■ u. 

a"* a a 

Now Eqs. fl45ll - fj48l) give closed differential equations for density and velocity 
perturbations, 5 and u. These are valid in the presence of the cosmological constant. 

In order to derive the axion contribution we have strictly ignored 
hH 


O * 

c^m 


order terms, and have kept up to only hrst order in 


27.2 






= 4.088 X 10"^^ 


10" W 


m 


IkpcV 

ao/k) 


(49) 


(50) 


In this section we have recovered c and h. The curvature perturbation variable 
ip is dimensionless with ip ~ where 5$ is the perturbed gravitational potential. 

In the Newtonian limit with c-goes-to-in£nity, the right-hand-sides of Eqs. fl4^ and 
fl47)l vanishes, thus recovering, except for the axion correction term, the well known 
Newtonian continuity and Euler equations valid to fully nonlinear order [19]; notice 
that in Einstein’s gravity these equations without ip corrections are valid only to the 
second order perturbation. 

The presence of Planck’s constant h in the axion correction term, despite of our 
classical analysis, can be traced to the presence of the constant in the potential 


VW = 2 


me 

h 


( 51 ) 


with a dimension [0] = (M/L)^/^. 
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The term involving m in Eq. (IT7|) is the pressure correction term arising from the 
axion nature of the system. The only difference of axion from the zero-pressure fluid 
(CDM) is represented by this axion pressure term. The left-hand-sides of Eq. fH6l) and 
(IT7|) . except for the axion pressure term, are identical to the Newtonian continuity and 
Euler equations, respectively [19]. The right-hand-side are third-order terms; these are 
pure general relativistic corrections in the zero-pressure fluid [T8| and are numerically 
negligible in the current paradigm of concordance cosmology for zero-pressure fluid [20] ; 
the same is true for axion in the super-Jeans scales, see below. 

The axion pressure contribution can be written as 


z'^^5Af = 


h^A AvT+J h^A Ay/g 


2 m?a‘^ y/l + 5 yfg ’ 

where Jl = gC = (1 -|- 6 )gC. Equations (15^ and (EZ]) show that the perturbed 

part of lapse function 6 Af is only loosely related to the perturbed pressure, see Eq. 
fl38p . We note that this so-called axion pressure term derived from general relativistic 
third order perturbation coincides (except for the scale factor) exactly with the one 
derived from Schrodinger equation in the non-relativistic limit in Minkowski space-time 
naE]. Although the authors of ^ emphasized that the result is due to Bose-Einstein 
condensate nature of the axion, we stress that our result is a purely classical one based on 
perturbation treatment of a massive scalar held in Einstein’s gravity without assuming 
non-relativistic limit; the resulting coincidence proves that the classical axion is non- 
relativistic indeed. 

Equations fHbll and fjTTll can be combined to give 


5 + 2-5- A'KGg5 + 4 [aV ■ (Ju)]' 
a 

= third order terms. 


■ (u ■ Vu) + 


h^A A^/m 

2iv?a‘^ vT+T 

(53) 


where third order terms in the right-hand-side are pure general relativistic terms all 
involving linear order ip term. Besides the axion pressure term, the left-hand-side is 
valid to fully nonlinear order in Newtonian theory (this is valid in the Einstein’s theory 
to the second order only). 

Comparing the gravity and pressure terms to the linear order we have the axion 
Jeans scale [21131 El [7| 


— 


2710 

kj 


Ggm?) 


5.4 X lO^^cmh 


m \ 

10-5evJ 


(54) 


with A = —Examination of Eq. fl53|) reveals that the same criteria applies to the 
third order. The axion pressure term has a role on scale smaller than the Solar System 
size for a canonical mass m ~ 10“®eV, thus completely negligible in all cosmological 
scales. This proves the CDM nature of the axion in the super-Jeans scale. 

Notice that the third order terms arising from the pure Einstein’s gravity in the 
right-hand-sides of Eqs. flT6l) . fITTI) and fl53|l include linear order ip term. From Eqs. 
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(USD-dlH]) we have 




^ a4c2m2a2'^' 

to the linear order. The homogeneous solution (which is constant in time) is the zero- 
pressure fluid part and the inhomogeneous solution is the additional contribution from 
the axion. From Eqs. (HHj) and (l55|) we have the fluid part and the additional axion part 
of ip as 


t\ 




fluid 




h^A 




(56) 


c2A ’ 4c2m2a2 

Thus, the pure axion part of p is negligible in the super-Jeans scale. We have 
‘^Ifluid = C'(x) = constant in time even in the presence of the cosmological constant, 
see Eq. fl55l) : p is the curvature perturbation in the comoving gauge which is the well 
known conserved quantity in the super-sound-horizon scale [22]. Even in sub-Jeans scale 
the pure axion part of p is simply decaying with (p|axion oc a~^ for A = 0. 

Therefore, in the cosmological scales we have proved that, to the third order 
perturbation, the axion can be identified as the CDM in the general relativistic context. 
In our previous works we have shown that the leading nonlinear power spectra (of 
density and velocity perturbations) in Einstein’s gravity (in the comoving gauge) give 
virtually identical results as in the Newtonian context izil. In the current paradigm 
of concordance cosmology the pure general relativistic contributions arising from the 
third order are entirely negligible in all scales compared with the relativistic/Newtonian 
results which are identified to the second order perturbation [10] . 

Using Eqs. flT3|) and flTSl) . Eq. fl26|l becomes 

{ (VA) u ■ VJ - ViA-V • [(VJ) u ■ VJ] }. (57) 


NO 

V- = 


ACm?a? 

This is the third-order vector-type perturbation generated from the axion. From Eqs. 
(jlj), (El) and (IT^ . to the third order, we have Tf /p = = Ui/a. Thus, we may identify 


Ui = ViU 


N) 


( 68 ) 


In case of the zero-pressure medium like the CDM, the nonlinear order vector-type 
velocity (rotational) perturbation is not generated from the scalar-type perturbation, 
and the homogeneous equation has decaying nature cx 1/a) in an expanding 

medium. This can be proved as the following. For a zero-pressure fluid in the comoving 
gauge, from the covariant momentum conservation in Eq. (3.9) of [12], we have 

(ah,^*^^)' -f = nonlinear terms all involving (59) 

valid to the fully nonlinear order. To the linear order, we have 6 J\f = 0 and cx 1/a. 
As the is a homogeneous and pure decaying mode we may set = 0; in 
fact, as we consider a vector-type perturbation generated from the pure scalar-type 
perturbation, we may ignore the linear order This can be continued to all higher 
order perturbations with SAf = 0 and cx 1/a; as the is homogeneous (not 
sourced by the scalar-type perturbation) and pure decaying mode, we may set = 0 
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to all orders in perturbation. [End of the proof.] Notice that in the matter-dominated 
era with = 0 = A, we have 5 oc a and u oc thus the axion-generated vector-type 
perturbation behaves as vf’'^ oc 

6. Discussion 

In this work we show that, in the super-Jeans scale, the axion behaves as the CDM 
to the third order perturbation. The axion pressure terms are important only in the 
cosmologically negligible scale for the canonical mass axion. In the axion-comoving 
gauge we have taken, in the super-Jeans scale, there exists the relativistic/Newtoninian 
correspondence to the second-order perturbations, and pure relativistic corrections 
appearing from the third-order are numerically negligible (in the current paradigm of 
concordance cosmology) j20]. Therefore, in those scales (including the super-horizon 
scale) the axion is indistinguishable from the CDM as a zero-pressure fluid. 

Concerning the axion pressure term, we show that our purely classical perturbation 
treatment of the axion as a massive scalar held in Einstein’s gravity gives an identical 
result from quantum mechanical treatment in the non-relativistic limit of the Bose- 
Einstein condensation m Although our result is valid only to third order 

perturbation, as the analysis is made in fully general relativistic context, the coincidence 
implies that the axion is non-relativistic indeed. In the regime where the axion pressure 
term becomes important, say for a extremely low mass axion where the axion Jeans scale 
becomes cosmologically relevant, besides the axion pressure term appearing from the 
linear order, the pure general relativistic contributions starting to appear from the third 
order through (p also have axion pressure contribution with decaying nature though, see 
Eq. (j56ll . 

We also show that the axion supports the vector-type velocity perturbation from the 
third order whereas the pressureless fluid without the vector-type velocity perturbation 
to the linear order does not generate the vector-type velocity perturbation to the 
nonlinear order. 

Our result includes the cosmological constant. Extending our analysis for a realistic 
cosmological situation in the presence of other components of fluids and helds is trivial. 
We emphasize that the spatial and temporal gauge condition in Eqs. (jb]) and fl^ and 
the special identihcation of perturbed velocity variable in Eq. (H3|l were essential to get 
the above equations. 

One missing contribution we ignored in this work is the simultaneously excited 
tensor-type (transverse-tracefree) perturbation. In our previous works [T8l |20] we have 
ignored the Y-term which comes from the vector-type perturbation generated by the 
nonlinear scalar-type perturbation. In the same sense the tensor-type perturbation is 
also generated by the nonlinear scalar-type perturbation. These scalar-generated vector- 
and tensor-type perturbations start to have their roles from the third-order perturbation, 
and their quantitative effects will be studied in a later work |21] . 

Whether the axion behaves as a CDM fluid to the fully nonlinear order in Einstein’s 
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gravity is left for future study; the fully nonlinear perturbation equations of the zero- 
pressure irrotational fluid in Einstein’s gravity are presented in Eqs. (56)-(60) in |llj . 
and our analysis up to Eq. (122|1 is valid to fully nonlinear order. 
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